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ABSTRACT 



We perform non-UTE calculations of lithium in late-type stars for a wide range of stellar parameters, including quantum mechanical 
cross-sections for collisions with neutral hydrogen and the negative hydrogen ion. Non-LTE abundance corrections for the lithium 
resonance line at 670.7 nm and the subordinate line at 610.3 nm, are calculated using ID MARCS model atmospheres spanning a 
grid Teff = [4000, 8000] K, \ogg = [1.0, 5.0], and [Fe/H] = [0.0, -3.0], for lithium'abundances in the range A(Li) = [-0.3,4.2]. 
The competing effects of ultraviolet over-ionization and photon losses in the resonance line govern the behaviour of the non-UTE 
effects with stellar parameters and lithium abundance. The size and sign of the non-LTE abundance corrections vary significantly over 
the grid for the 670.7 nm line, but are typically positive and below 0.15 dex for the 610.3 nm, line. The new collisional data play a 
significant role in determining the abundance corrections. 
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1. Introduction 

Stellar lithium abundances are continuing to attract strong 
interest among the astronomical community. In particular, 
many studies are devoted to constraining and finding expla- 
nations for the famous "Spite plat eau" of warm m etal-poor 
halo stars, first discovered by Spi te & Spitd (|1982[). For re- 
cent s tu dies of lithium in hal o field st ars see e.g. Ryan et al] 



(1200 lb. ICharbonnel & Priiiiasl (l2005h. lAsplund et all (l2006h. 



iBonifacio etall ( l2007h . and iHosford et alj (l2009l) . To put con- 
strains! on the primordial lithium abundance in the Universe, 
it is of fundamental importance to determine if and by how 
much the Li abundances increase with increasing metallicity on 
the Spite plateau. The behaviour of Li abundance with effec- 
tive temperature is of equal importance, as it could lend sup- 
port to the notion that atomic diffusion is acting in metal-poor 
stellar atmospheres. Lithium depletion through atomic diffusion 
has been suggested as a solution to the discrepancy between 
the Spite plateau abundance and th e predicted value of the pri- 
mor dial lithium abundance (see e.g. lKorn et al.ll20O7l : lLind et alj 
12009 ). To accurately infer both the mean lithium abundance 
in the halo and the abundance behaviour with metallicity and 
effective temperature, it is crucial to have a realistic descrip- 
tion of the lithium line formation. Previous non-LTE an alyses 
spanning a large stellar-parameter space ( Carlsson et al| 1994t 
iPavlenko & Magazzu 1996; Takeda & Kawanomotai2005 l) have 
shown that departures from LTE are generally small but signifi- 
ca nt at the required accu racy. 

iBarklem et al.l ( l2003i) applied quantum-mechanical calcula- 
tions of cross sections for inelastic collisions with neutral hy- 
drogen (iBelvaev & BarklemI 120031: ICroft et alj [1999) in non- 
LTE analysis of the Sun and two metal-poor stars. Their find- 



ings point to a negligible influence on the statistical equilib- 
rium of lithium from collisional bound-bound transitions with 
hydrogen, but a significant influence from charge transfer reac- 
tions, specifically mutual neutralization and ion-pair production 
(Li* -H H t5 Li^ + H ). Other non-LTE investigations have re- 
lied on the classical Drawin re cipes (jDrawin 1968), as given by 
ISteenbock & Holwe geJ (119 84 and iLambert (1993), or the free 
electron model of rKaulakv3 Tl985l) . for estimates of collisions 
with neutral hydrogen. None except Barklem et al. have included 

the influential charge transfer reaction.^ 

We have extended the study by iBarklem et al.l (|2003|) to 
cover a large cool-star grid and calculate non-LTE abundance 
corrections for the lithium 670.7 nm (2s-2p) and 610.3 nm (2p- 
3d) lines. 

2. Setup 



* Tables 2 and 3 are only available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
|http://cdsweb.u-strasbg.fr/ , 



We use the r adiative-transfer code MULTI, version 2.3 dCarlssonl 
1198611199 2') to perform non-LTE calculations. The model atom 
used include s the same 20 e nergy levels for neutral lithium as 
described in ICarlsson et al.l (1 19941) . plus the Li II ground state. 
The highest considered level in Li I has principle quantum num- 
ber n - 9. We have used TOPbase data (Peach et al. 1988) for 
energy levels, oscillator strengths, and photo-ionization cross- 
sections for levels with orbital quantum number I < 3. For 
the remaining levels, hydrogenic values are used. For the res- 
onance line at 670.7 nm , we ado pt the oscillator strength / = 
0.7468, as calculated bv I Yan etal. (1998), and consider six hy- 
perfine componen ts in ^Li (neglecting ^Li) , with measured wave- 
lengths given by iSansonetti et al.l (Il995h . For the subordinate 
line at 610.3 nm, we adopt / = 0.6386, also determined by 
Yan et al, and account for the three fine-structure components 
with wavelengths determined bv .Lind gard & Nielson (1977|. 
For both lines, van-der -Wa als-broadening parameters follow 
lAnstee&O'JVIaral (119951) and lBarklem & 0']VIaral(ll997l) . Stark 
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broadening is unimportant in the late-type atmospheres of inter- 
est here and is therefore neglected. Table 1 lists wavelenghts, os- 
cillator strengths, and broadening data adopted for the 670.7 nm 
and 610.3 nm lines. 

Cross-sections for collisional excitation by electrons are 
taken fromJPark (1971) and collisiona l iori i zation by elec- 
trons from Seatoni (119621) as given by AUenl (11976 '). We add 
rate coefficients for excitation and de-excit ation from collisions 
with n eutral hydrogen atoms acco rding to iBelvaev & BarklemI 
(l2003h and iBarklem et all (l2003h . as well as charge transfer 
reactions with n eutral hydrogen an d the negative hydrogen 
ion according to ICroft et al] (Il999h . The number abundance 
of the negative hydrogen ion is calculated by assuming LTE. 
Ionization by hydrogen atom impact is not included. For low- 
lying states at low collision energies, it is expected to be negli- 
gibl e compared to charge t ransfer reactions and excitation (see 
e.g. lKrstic & Schultz 2009). 

A grid of ID, LTE, opacity-s ampling, MARCS model at- 
mospheres (IGustafsson et al.ll2008h is used in the analysis. Non- 
LTE computations for lithium as a trace element are performed 
in the plane-parallel approximation, for models with Teff - 
[4000, 8000] K, log ^ = [1.0,5.0], [Fe/H] = [0.0, -3.0], for 
lithium abundances in the range A(Li) = [-0.30,-4.20]. The 
highest effective temperature is 5500 K for models with log g - 
1.0, 6500 K for log g = 2.0, 7500 K for log g = 3.0, and 8000 K 
for \ogg > 4.0. For models with logg > 3.0, we adopted a mi- 
croturbulence parameter - l.Okms"' and - 2.0kms"', 
and for models with logg < 3.0, we adopted = 2.0km s ' and 
ft = 5.0km s"'. In total, 392 atmospheric models are included in 
the computations. 

In the statistical equilibrium calculations, the background 
opacities include continuum and line opacities provided by the 
MARCS model atmospheres. MULTI thus calculates a line- 
blanketed photo-ionizing radiation field while solving for the 
statistical equilibrium of lithium. Based on the LTE and non- 
LTE equivalent widths obtained for the stellar-parameter grid, 
we subsequently define a non-LTE correction for each abun- 
dance point as the difference between its LTE lithium abundance 
and the non-LTE abundance that corresponds to the same equiva- 
lent width. The limiting equivalent widths for which corrections 
are given are set to 0.01 pm and 100pm. The equivalent width is 
obtained by numerical integration over the line profile, consider- 
ing a spectral region that extends ±9 nm from the line centre. For 
lines weaker than 50 pm the numerical precision is better than 
0.01 pm whereas the equivalent widths of the strongest lines are 
determined to within 0. 1 pm. 

3. Results 

The abundance corrections and LTE and non-LTE equivalent 
widths for each grid point are given in Table 2 for the 670.7 nm 
line, and in Table 3 for the 610.3 nm line. In the following we 
summarize the most important non-LTE effects and describe the 
size of the corrections over the cool- star grid. Because ou r re- 
sults are in qualitative agreement with lCarlsson et all d 19941) . we 
refer the reader to this paper for a detailed explanation of the 
non-LTE line formation of lithium. 

At A(Li) < 2 for the coolest giants and A(Li) < 3 for the 
hottest giants, the dominant non-LTE effect is over-ionization 
of neutral lithium, driven mainly by the super-thermal radiation 
field at A ^ 349.5 nm (Jy > By, where Jy is the mean inten- 
sity and By the Planck function), corresponding to the photo- 
ionization edge for the first excited level. This leads to smaller 
number populations of the ground (2s) and first excited (2p) 



Table 1. Wavelengths, oscillator strengths, and broadening data 
for the two considered lines. 
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r [rad 5-'] is the natural broadening parameter. 

<T [a.u. ] is the broadening cross-section for collisions 
with neutral hydrogen at relative velocity v = 10''m s-'. 
lAnstee & O'Mara 1995) 
''^^ a is the velocity dependence of cr. 



states compared to LTE and in higher values of the line source 
functions. Both the effect on the source functions and the loss of 
line-opacity weaken the lines in non-LTE, making the abundance 
coiTections positive. The coiTections range from being close to 
zero for the hottest giants to a maximum -1-0.45 dex for cool, 
metal-rich giants for the 670.7 nm line. For the 610.3 nm line, 
the corrections are slightly smaller, < -1-0.3 dex. 

In dwarfs, the Jy - By excess is smaller in the 2p photo- 
ionization continuum and the radiation field is sub-thermal, 
Jy < By, in other influential continua, e.g. those correspond- 
ing to photo-ionization from 3p and 3d. In layers close to con- 
tinuum optical depth unity, hot dwarfs have overpopulations in 
the ground state, which in combination with a decrease in the 
resonance-line source function strengthens the 670.7 nm line. 
Abundance corrections are therefore negative or approximately 
zero for this line for dwarfs hotter than 6000 K. 

At A(Li) > 2 for the coolest stars and A(Li) > 3 for the 
hottest stars, photon losses in the resonance line become appar- 
ent, driving recombination from Li II in both dwarfs and giants. 
This mainly affects the population of the ground state and to 
a lesser degree higher excited states. The resonance-line source 
function drops accordingly and the non-LTE corrections become 
less positive or more negative with increasing abundance over 
the whole grid. For metal-poor, extremely lithium-rich giants, 
the corrections can reach almost -1.0 dex. 

Since the number population of the first excited state and the 
source function of the subordinate line are not as severely af- 
fected by the photon losses in the resonance line, the 610.3 nm 
abundance coiTections for A(Li) < 4.0 are approximately con- 
stant for fixed stellar param eters. 

As described in Barklet iTet al.l (l2003l) . including hydrogen 
collisions strengthens the collisional coupling between singly 
ionized lithium and neutral lithium, especially through the 
charge transfer reaction involving H and Li I in the 3s state 
and its inverse reaction. The number populations of low excited 
states increase when charge transfer is included, while the line 
source functions are barely affected. In turn, the abundance cor- 
rections become generally lower, i.e. less positive or more nega- 
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Fig. 1. Non-LTE abundance corrections for 670.7 nm as func- 
tions of LTE lithium abundance for the indicated stellar pa- 
rameters in the lower left corner of each plot, to be read as 
Teff/ logg/[Fe/H]. Solid: Including charge transfer reactions 
with hydrogen as well as bound-bound transitions due to col- 
lisions with neutral hydrogen. The cross-sections are calculated 
with quantum mechanics (Belyaev & Barklem 2003; Croft et al. 
1999). Dotted: Neglecting charge transfer reactions, but includ- 
ing bound-bound transitions, calculated with quantum mechan- 
ics. Dashed: Neglecting charge transfer reactions, but includ- 
ing hydrogen collisions for bound-bound and bound-free transi- 
tions. The cross sections are calculated with the classical Drawin 
recipe. 

tive. Notice that the added colUsions do not always have a ther- 
malising effect, but can contribute to over-population of levels 
compared to LTE. The effect is illustrated in Fig. 1, for a selected 
metal-poor and solar-metallicity dwarf and giant. As seen in Fig. 
1, the change in A(Li)non-LTE -^(Li)LTE when including charge 
transfer is -0.12 dex for the metal -poor giant and -0.06 dex for 
the metal-poor dwarf. 

While the charge transfer reaction has a significant influence 
on the statistical equilibrium of lithium, including bound-bound 
(and bound-free) transitions due to collisions with neutral hydro- 
gen has no influence at all. However, this is not true if one relies 
on the classical Drawin recipe for estimates of the cross-section 
for those collisions. For comparison. Fig. 1 also shows the results 
obtained when neglecting charge transfer reactions but includ- 
ing excitation and ionization by neutral hydrogen according to 
Drawin's recipe, with the rate formula given by Lambert C 1 993i) 
and with the scaling factor Sh - 10. Generally, the Drawin 
cross-sections are much higher than the quantum mechanical 
calculations by Belyaev & Barklem. When adopting the higher 
classical rates, the ground state is more populated and the abun- 
dance corrections become lower, by 0.05 dex, compared to 

the results obtained with quantum mechanical rates. 



Fig. 2. Panels as in Fig. 1. Solid: Same as in F ig. 1 . MARCS 
model atmospheres from lGustafsson et al.l (12008 ^ are used, with 
= 2kms~'. Dashe d: MARCS model atmospheres from 
iGustafsson et al.l (1 19751) are used, with — 2kms ~ '. Do tted: 
MARCS model atmospheres from IGustafsson et al.l (l2008l) are 
used, with microturbulence parameter set to = 5 km s^' in the 
two upper panels and = 1 km s ' in the lower panels. 



Table 2. Non-LTE abundance corrections and equivalent widths 
for the Li I /i = 670.7 nm line. The whole table can be retrieved 
in machine-readable format from the CDS. 
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Table 3. Non-LTE abundance corrections and equivalent widths 
for the Li I /i = 610.3 nm line. The whole table can be retrieved 
in machine-readable format from the CDS. 
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4. Discussion 

We h ave compared our results for the lithium resonance line to 
ICar lsson et al. ( 1994) abundance corrections for dwarfs and gi- 
ants with [Fe/H] = [-3.0,0.0] to iPavlenko & Magazzul d 19961) 
results for dwarfs and subgiants with solar metallicity and to 



iTakeda & Kawanomotd (l2005h corrections for dwarfs and sub- 
giants with [Fe/H] = [0.0, -1.0] (see Fig. 3). To aid in the com- 
parison we converted the LTE and non-LTE equivalent widths 
listed in Table 1 in Paylenko & Magazzu (1996) to abundance 
corrections. The l argest difference is 0.1 2 dex between our re- 
sults and those of iPavlenko & Magazzul and 0.15 dex between 
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Fig. 3. Panels as in Fig. 1 . S olid: Same as in Fig. 1 
and 2. Dotted: Carlsson et al.l (11994') corrections. Dashed: 
Takeda & Kawanomotol (l2005l) corrections. Dash - dotted: 
Pavlenko & Magazzul d 19961) corrections. 



our results and those of iTakeda & Kawanomotol but ge nerally 
the values agree quite well. The ICarlsson et al.l d 19941) values 
agree with ours to within 0.15 dex for dwarfs and 0.20 dex 
for giants. Such differences are reasonable considering dif- 
ferences in th e mo del atom and in the atmospheric models. 
ICarls son et"aD use 'Gustaf sson et all (11975,) MARCS models 
and Takeda & Kawanomotc^ (l2005h and iPavlenko & Magazzul 
([1996) use K uruczi (|1993) ATLAS9 models. The effects on the 
non- LTE abundance c orrections when using MARCS mo dels 
from iGustafsson et al.1 d 19751) and iGustafsson et al.1 d2008l) are 
illustrated in Fig 2. Especially for the two giants, the choice of 
model atmosphere is important for the corrections. The newer 
models include more line opacity, causing a steeper temperature 
gradient in the upper part of the photosphere, which increases 
the ultraviolet Jy - By excess and leads to more over-ionization. 
The use of newer model atmospheres thus leads to more positive 
or less negative abundance corrections and partly cancels with 
the effect of including hydrogen collisions (see Sect. 3) 

Varying the microturbulence between 1 km s"' and 2 km s"' 
barely affects the abundance corrections (lower panels in Fig. 2), 
especially for low lithium abundances. In a giant star, the non- 
LTE abundance corrections become systematically higher when 
adopting a microturbulence of 5 kms"' instead of 2kms"', the 
differences being especially significant when the line is strong 
(upper panels Fig. 2). However, the parameter still has some 
significance at low lithium abundances, when the formation of 
the line itself is unaffected by the microturbulence. This is be- 
cause the choice of microturbulence influences the amount of 
line opacity included in the computations of the model atmo- 
sphere and consequently also the atmospheric temperature gra- 
dient, which in turn partly drives the departures from LTE. 

Our results are valid within the assumptions of ID model 
atmospheres in LTE and hydrostatic equilibrium. Non-LTE cal- 
culations for lithium in 3D, hydrodynamical model atmospheres 
have to this date bee n performed only for a few types of stars 
dAsplund et al.ll2003l) . and we plan to extend our work to 3D in 
the future. 
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